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Antioxidant  capacity  of  3,4-dihydroxybenzaldehyde 
and  the  hydroperoxyl  radical  in  solution: 
Atheoretical  approach 


Abstract 

Herein,  we  studied  the  antioxidant  capacity  of  the  3,4-dihydroxybenzaldehyde  against 
the  hydroperoxyl  radical  by  means  of  a  thermochemical  and  kinetic  analysis  using  the 
Theory  of  Functional  Density.  The  energies  of  the  reaction  barriers  were  computed  at 
the  M05-2X/6-31 1  G(d,p)  level  of  theory  employing  the  Transition  State  Theory.  The 
thermochemical  analysis  identified  the  hydrogen  transfer  mechanism  as  the  main  route 
taken  by  the  reaction.  The  kinetic  analysis  indicates  that  3,4-dihydroxybenzaldehyde  is 
effective  to  inhibit  this  radical  since  the  overall  kinetic  constant  was  higher  than  the 
oxidation  constant  of  fatty  acids  present  in  cell  membranes. 
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Resumen 

Esta  investigación  estudió  la  capacidad  antioxidante  del  3,4-dihidroxibenzaldehído, 
frente  al  radical  hidroperoxilo  por  medio  de  un  análisis  termoquímico  y  cinético: 
se  utilizó  la  Teoría  de  los  Funcionales  de  la  Densidad.  Las  energías  de  las  barreras  de 
las  reacciones  se  calcularon  recurriendo  a  la  Teoría  del  Estado  de  Transición,  con  el 
nivel  de  teoría  M05-2X/6-31 1+G(d,p).  Los  resultados  termoquímicos  muestran  que  el 
mecanismo  de  transferencia  de  hidrógeno  es  la  vía  principal  que  toma  la  reacción.  Los 
cálculos  cinéticos  evidencian  que  el  3,4-dihidroxibenzaldehído  es  efectivo  para  inhibir 
a  este  radical,  ya  que  la  constante  de  velocidad  global  resultó  mayor  a  la  constante  de 
oxidación  de  los  ácidos  grasos  presentes  en  las  membranas  celulares. 

Palabras  clave:  Cinética  Química,  estrés  oxidativo,  radical  libre,  capacidad  antioxidante. 
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INTRODUCCIÓN 

La  Cinética  Química  se  encarga  del  estudio  de  la  velocidad  con  la  que  ocurren  las 
reacciones  químicas  y  posee  una  relación  intrínseca  con  el  mecanismo  que  las  genera. 
Los  modelos  cinéticos  permiten  la  comprensión  de  procesos  químicos  importantes 
como  la  descomposición  de  los  alimentos,  el  crecimiento  de  microorganismos,  la 
disgregación  del  ozono  estratosférico  o,  en  este  caso,  la  determinación  de  la  capacidad 
antioxidante  de  un  compuesto  [1,  2].  Desde  el  punto  de  vista  teórico,  los  cálculos 
cinéticos  se  realizan  al  conocer  las  funciones  de  partición  para  las  sustancias  que 
participan  en  la  reacción;  entonces  es  posible  calcular  la  constante  de  velocidad,  si  se 
emplea  la  teoría  del  estado  de  transición  (TST,  por  sus  siglas  en  inglés)  [3, 4].  El  método 
para  determinar  la  capacidad  antioxidante  de  un  compuesto  consiste  en  comparar 
la  constante  de  velocidad  del  proceso,  con  la  velocidad  de  oxidación  de  los  lípidos 
presentes  en  las  membranas  celulares.  La  oxidación  lipídica  se  genera  por  las  especies 
reactivas  de  oxígeno  durante  el  estrés  oxidativo. 

1.  Estrés  oxidativo 

El  estrés  oxidativo  (OS)  es  un  proceso  biológico  que  se  desarrolla  en  disolución  y  ocurre 
en  los  organismos  vivos  como  consecuencia  del  metabolismo  aerobio.  El  OS  da  lugar  a 
un  entorno  que  se  produce  por  el  desequilibrio  entre  la  producción  de  especies  reactivas 
de  oxígeno  (ROS)  y  la  capacidad  del  sistema  biológico  de  reparar  el  daño  inducido  por 
estas  especies.  Parte  de  los  ROS  se  generan  por  procesos  fisiológicos  que  se  derivan 
del  oxígeno  molecular  [5],  el  cual  produce  especies  como  el  radical  anión  superóxido 
(02-),  el  radical  hidroxilo  (HO-),  el  oxígeno  singulete  (]02),  radicales  alcoxilo  (RO),  radicales 
peroxilo  (ROO-)  y,  además,  genera  especies  reactivas  de  nitrógeno  [6].  En  las  células,  la 
principal  fuente  de  producción  de  ROS  es  la  cadena  respiratoria  mitocondrial  que  utiliza 
aproximadamente  el  80-90%  del  02  que  consume  una  persona  [7].  La  producción  de  ROS 
también  se  da  por  la  exposición  a  la  radiación  ionizante,  radiación  UV  o  metabolismo  de 
algunos  fármacos  [8]. 

Los  ROS  regulan  varios  procesos  celulares  en  los  mamíferos,  como  la  producción  de 
hormonas  de  crecimiento,  la  comunicación  entre  células,  la  unión  de  las  proteínas  G,  entre 
otros  [9].  Cuando  el  nivel  de  especies  reactivas  de  oxígeno  supera  significativamente  la 
capacidad  de  las  defensas  antioxidantes  en  el  sistema  biológico,  estas  tienen  efectos 
tóxicos  en  el  medio  [10].  La  sobreproducción  de  ROS  contribuye  a  la  patogénesis  de 
numerosas  enfermedades  neurodegenerativas  relacionadas  con  la  edad  y  el  deterioro 
cognitivo  en  los  seres  humanos,  debido  a  su  ataque  en  los  ácidos  grasos  insaturados  de 
las  membranas  celulares  y  las  cadenas  de  ADN  [1 1]. 

1.1.  Reactividad  de  las  especies  de  oxígeno 

La  reactividad  de  los  ROS  se  evidencia  por  medio  del  potencial  de  reducción;  se  conoce 
que  cuanto  mayor  es  el  potencial  de  reducción  de  los  ROS,  mayor  es  su  capacidad  de 
oxidación.  El  peróxido  de  hidrógeno  no  se  involucra  directamente  en  el  inicio  de  la 
oxidación  de  los  lípidos,  debido  a  que  su  potencial  de  reducción  de  320  mV  es  inferior  a 
los  600  mV  de  los  ácidos  grasos  insaturados.  El  potencial  de  reducción  del  radical  anión 
superóxido  es  de  940  mV:  no  es  lo  suficientemente  fuerte  como  para  extraer  hidrógeno 
de  los  ácidos  grasos  insaturados  [1 2].  El  radical  hidroxilo  tiene  un  potencial  de  reducción 
estándar  más  elevado  (2310  mV),  pero  este  se  genera  a  partir  del  H202  [13].  El  radical 
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hidroperoxilo,  con  1060  mV,  es  más  reactivo  que  el  radical  anión  superóxido  [14]  y,  por 
ello,  su  presencia  en  la  membrana  mitocondrial  de  las  células  origina  la  peroxidación 
de  lípidos.  La  tasa  de  extracción  de  H-  de  los  ácidos  grasos  insaturados  varía  entre  1-3  x 

1 03  M_1  s~\  lo  que  representa  un  paso  crítico  en  el  proceso  de  peroxidación  lipídica  [1 2]. 

7. 7. 7  Formación  del  radical  hidroperoxilo 

El  radical  anión  superóxido  se  forma  en  la  cadena  de  transporte  de  electrones  de  la 
membrana  mitocondrial,  a  partir  del  oxígeno  triplete.  En  la  célula,  cuando  el  anión 
superóxido  formado  se  difunde  hacia  el  medio  acuoso  celular,  forma  su  ácido  conjugado: 
el  radical  hidroperoxilo,  según  la  reacción  (1)  [15]. 

02*-(ac)  +  H+(ac)  HOO*(ac)  (1) 

Esta  reacción  muestra  que  los  dos  radicales  se  encuentran  en  equilibrio  en  medio 
acuoso.  La  reacción  inversa,  correspondiente  a  la  deprotonación  del  HOO-  tiene  una 
constante  de  equilibrio,  «H00  =  1,6  x  10~5  mol/L  [16].  Esto  indica  que,  en  medio  acuoso,  la 
formación  del  HOO-  es  favorable. 

El  radical  hidroperoxilo  se  difunde  a  través  de  la  membrana  hidrofóbica  [1 1]  y  produce 
una  variedad  de  radicales  lipidíeos  (R-)  debido  a  su  capacidad  de  extraer  hidrógeno 
de  los  ácidos  grasos  presentes  en  las  células,  lo  que  da  lugar  a  la  aparición  de  varias 
enfermedades  [17,  18].  Las  velocidades  de  reacción  del  radical  hidroperoxilo  con 
los  ácidos  linoleico,  linolénico  y  araquidónico  son  1,2  x  103,  1,7  x  103  y  3,1  x  103  IVL1 
s_1,  respectivamente  [12,  19].  Estas  constantes  son  relativamente  altas  y,  por  tanto, 
determinantes  en  los  procesos  oxidativos. 

Algunas  sustancias  químicas,  como  los  compuestos  fenólicos,  tienen  la  capacidad  de 
inhibir  el  daño  causado  por  estos  radicales  y  ejercer  un  efecto  protector  a  nivel  celular  [20]. 

2.  Actividad  antioxidante  de  los  compuestos  fenólicos 

Los  compuestos  fenólicos  son  una  clase  de  antioxidantes  que  constan  de  uno  o 
más  grupos  hidroxilo  (HO-)  unidos  directamente  a  uno  o  más  anillos  aromáticos.  Los 
polifenoles  se  distribuyen  ampliamente  en  las  plantas,  frutas,  verduras,  el  vino  y  el  té 
[21, 22].  Los  compuestos  fenólicos  poseen  la  propiedad  de  prevenir,  reducir  o  reparar  el 
daño  inducido  por  las  especies  radicales  de  oxígeno  y  nitrógeno  en  las  biomoléculas. 

El  efecto  inhibidor  de  los  antioxidantes  polifenólicos  sobre  los  ROS  se  desarrolla  por 
medio  de  mecanismos  que  surgen  de  la  reacción  directa  con  radicales  libres  [23]. 
Estos  son  la  formación  radicalaria  de  aducios  (RAF),  transferencia  de  hidrógeno  (HT)  y 
transferencia  de  electrón  (SET). 

En  el  mecanismo  de  transferencia  de  hidrógeno,  los  antioxidantes  inactivan  los  radicales 
libres  por  medio  de  la  transmisión  del  átomo  de  hidrógeno.  Así,  un  antioxidante 
polifenólico  genérico,  ArOH,  reacciona  con  el  radical  libre  R-  y  transfiere  un  átomo  de 
hidrógeno  a  través  de  la  ruptura  homolítica  del  enlace  O-H,  para  generar  radicales  más 
estables,  como  se  representa  en  la  reacción  (2). 
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Los  productos  de  la  reacción  conducen  a  la  formación  de  las  especies  RH  inofensivas  y 
el  radical  oxidado  ArO-.  El  radical  que  se  forma  es  menos  reactivo  con  respecto  a  R-  [24]. 
El  mecanismo  de  transferencia  de  electrón  permite  la  donación  de  un  electrón  por  parte 
del  antioxidante  ArOH  al  radical  R-  según  la  reacción  (3). 


Ar  Ar 


El  anión  R-  es  una  especie  con  un  número  par  de  electrones,  mientras  que  el  radical 
catlónico  ArOH+-  es  una  especie  radical  menos  reactiva.  En  particular,  el  ArO-  y  ArOH-+ 
son  estructuras  aromáticas  en  las  que  el  electrón  impar  -originado  por  las  reacciones 
con  el  radical  libre-  tiene  la  posibilidad  de  extenderse  por  toda  la  molécula  y  producir 
una  estabilización  radical  por  resonancia  [23]. 

Otro  mecanismo  que  puede  ocurrir  entre  un  antioxidante  y  la  especie  radical  R-,  es  la 
formación  de  aductos  de  acuerdo  con  la  reacción  (4).  Esta  es  importante  en  antioxidantes 
que  poseen  dobles  enlaces  en  su  estructura,  como  en  el  caso  de  los  carotenoides  [25]. 


Ar  Ar 


El  3,4-dihidroxibenzaldehído  (nombre  IUPAC)  o  aldehido  protocatéquico  (PAL)  es  un 
compuesto  de  tipo  fenólico  (Fig.  1)  estructuralmente  similar  a  otras  moléculas  con 
capacidad  antioxidante  [26].  Estudios  experimentales  revelaron  que  este  compuesto 
inhibe  la  producción  de  especies  reactivas  de  oxígeno  relacionadas  con  la  aparición  de 
la  enfermedad  de  Parkinson  [27]. 
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FIGURA  1 .  Estructura  y  numeración  de  los  átomos  del  3,4-dihidroxibenzaldehído.  El  3,4-dihidroxibenzaldehído  está 
presente  en  la  salvia miltiorrhiza  [28]  (una  de  las  plantas  de  la  medicina  tradicional  China)  y  otras  plantas  y  semillas 
comunes  como  la  cebada.  Jeong  etal.  purificaron  e  identificaron  este  compuesto  en  la  cebada  y  evaluaron  sus  propie¬ 
dades  quimiopreventivas  contra  el  cáncer,  a  través  de  la  inhibición  del  daño  oxidativo  del  ADN  [29].ativo  del  ADN  [29]. 


2.1.  Cuantificación  de  la  actividad  antioxidante 

La  actividad  antioxidante  puede  expresarse  como  equivalentes  de  trolox,  una  unidad 
de  medida  relativa  con  respecto  al  trolox  (un  antioxidante  de  referencia)  [30]  o  como 
porcentaje  de  radicales  libres  removidos  o  inhibidos  [29].  En  un  estudio  experimental  o 
teórico  se  determina  la  capacidad  antioxidante  al  saber  que,  si  la  constante  de  velocidad 
de  la  reacción  de  inhibición  de  radicales  es  mayor  a  la  velocidad  de  extracción  de  H-  de 
los  ácidos  grasos  insaturados  por  el  radical  hidroperoxilo,  el  compuesto  es  efectivo  para 
la  eliminación  o  inhibición  de  estas  especies  [2, 31]. 


MATERIALES  Y  MÉTODOS 

Los  cálculos  electrónicos  de  las  especies  químicas  se  realizaron  empleando  la  Teoría  de 
los  Funcionales  de  la  Densidad  (DFT).  En  particular,  se  utilizó  el  funcional  híbrido  M05-2X 
[32],  diseñado  para  cálculos  computacionales  en  termoquímica  y  cinética  [33].  Según 
Truhlar  et  al.,  aquel  se  encuentra  entre  los  más  efectivos  para  calcular  las  energías  de 
reacción  que  involucran  radicales  libres  [32];  al  considerar  que  estos  sistemas  son  de 
capa  abierta,  [34]  se  empleó  el  funcional  correspondiente  UM05-2X. 

El  conjunto  de  funciones  base  que  se  manejó  es  el  6-31 1+G(d,p)  [35]:  un  conjunto 
de  bases  de  valencia  dividida  con  funciones  de  polarización  (d,p)  para  la  descripción 
del  estado  basal  de  cada  átomo  [36].  El  signo  +  indica  que  la  base  se  complementa 
con  funciones  difusas  en  los  átomos  pesados,  para  permitir  a  los  orbitales  ocupar  un 
mayor  espacio  y  describir  mejor  a  los  electrones  que  se  encuentran  distantes  del  núcleo. 
Estas  funciones  son  básicas  para  representar  apropiadamente  a  especies  ionizadas  y 
moléculas  con  pares  de  electrones  no  compartidos. 
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Para  incluir  los  efectos  del  disolvente,  se  empleó  el  modelo  de  disolvente  universal, 
SMD  [37]  que  permite  el  cálculo  de  energías  libres  de  solvatación  para  solutos  neutros  e 
iónicos  en  agua  y  en  disolventes  no  acuosos.  En  la  presente  investigación,  se  realizaron 
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cálculos  con  dos  disolventes:  agua  y  etanoato  de  pentilo,  para  imitar  el  medio  acuoso  y 
no  polar  [38, 39]  de  las  células.  El  agua  presenta  una  constante  dieléctrica,  e=  78,3553  y 
el  etanoato  de  pentilo  £=4,7297. 

Los  mínimos  locales  y  los  estados  de  transición  se  identificaron  por  el  número  de 
frecuencias  imaginarias,  0  y  1,  respectivamente.  Para  verificar  que  los  estados  de 
transición  encontrados  son  los  correctos,  se  efectuó  el  cálculo  de  Coordenada  Intrínseca 
de  Reacción  (IRC)  [40].  Este  es  un  cálculo  de  trayectoria  de  mínima  energía  de  reacción 
que  permite  verificar  si  el  estado  de  transición  se  conecta  con  los  reactivos  y  productos 
propuestos  [41]. Todos  estos  cálculos  se  efectuaron  con  Goussion  09  [42]. 

Para  determinar  la  capacidad  antioxidante  por  medio  de  la  Cinética  Química  se  empleó  el 
programa  Eyringpy,  un  programa  escrito  en  Python  específico  para  calcular  propiedades 
termoquímicas  y  cinéticas  [43].  Aquí  se  evaluaron  los  caminos  de  reacción  HT,  RAF  y  SET 
que  tomaría  el  3,4-dihidroxibenzaldehído  frente  al  HOO. 


Las  constantes  de  velocidad  se  calcularon  de  acuerdo  con  la  teoría  del  estado  de 
transición  para  los  mecanismos  RAF  y  HT.  LaTST  es  una  formulación  mecano-estadística 
para  computar  las  constantes  de  velocidad  de  las  reacciones  químicas  [3, 4];  se  basa  en  el 
supuesto  de  que  una  reacción  química  avanza  a  lo  largo  de  una  ruta  de  mínima  energía 
que  conecta  reactivos  y  productos  (Fig.  2).  La  configuración  que  divide  la  superficie  de 
energía  potencial  (PES)  en  regiones  de  reactivos  y  productos,  se  conoce  como  estado  de 
transición  (TS)  y  es  un  máximo  en  la  dirección  de  la  coordenada  de  reacción  [44].  La  TST 
supone  que  existe  un  quasi-equilibrio  entre  los  reactivos  y  el  estado  de  transición  [45]  de 
donde  se  deriva  la  Ec.  (5), 


Je  — 

^TST 


kBT  -ag^'m,rt 


(5) 


en  la  cual  kB  es  la  constante  de  Boltzmann,  T  es  la  temperatura,  h  es  la  constante  de 
Planck  y  R  es  la  constante  universal  de  los  gases.  Las  unidades  de  k^  son  L  mol-1  sT 


FIGURA  2.  Representación  unidimensional  de  la  energía  de  Gibbs  de  una  reacción. 
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En  las  reacciones  que  involucran  transferencia  de  hidrógeno  puede  ocurrir  el  efecto 
túnel  [45]  Su  importancia  en  Cinética  Química  radica  en  que  la  velocidad  de  este  tipo 
de  reacciones  podría  aumentar  varios  órdenes  de  magnitud  [46].  El  factor  que  permite 
conocersiexistetunelamientoen  una  reacción  química  eselcoeficientedetunelamiento 
o  coeficiente  de  transmisión,  el  cual  se  introduce  como  un  factor  multiplicativo,  k,  en  la 
Ec.  (5).  En  esta  investigación,  algunas  de  las  constantes  de  velocidad  se  corrigieron  por 
efecto  túnel  mediante  la  aproximación  de  Eckart  [47],  si  k  resulta  mayor  que  1  existe 
tunelamiento  y  si  resulta  menor  que  1 ,  se  lo  descarta.  La  ecuación  que  incluye  tal  factor 
se  representa  en  la  Ec.  (6). 

kTsT=Kb¡Le-^  (6) 

h 

En  el  caso  de  las  reacciones  de  transferencia  de  electrón  no  ocurre  la  formación  de  un 
estado  de  transición,  debido  a  que  las  moléculas  de  disolvente  no  están  en  equilibrio 
con  las  cargas  de  las  especies  que  se  oxidan  y  se  reducen  [48].  En  este  tipo  de  reacciones, 
las  moléculas  de  disolvente  se  reorganizan,  es  decir  se  alejan  o  acercan  de  acuerdo  al 
nuevo  tamaño  de  las  moléculas  (o  iones). 

Marcus  propuso  una  configuración  en  la  que  las  energías  de  Gibbs  de  los  reactivos  y 
productos  son  iguales  [48, 49];  así,  la  expresión  matemática  que  permite  el  cálculo  de  la 
constante  de  velocidad  para  las  reacciones  en  las  que  ocurre  transferencia  electrónica, 
se  expresa  por  medio  de  la  Ec.  (7), 

kET=^Le-^r'RT  (7) 

h 

donde  AG*ET  es  la  energía  de  Gibbs  de  activación. 

Además,  las  reacciones  de  transferencia  de  electrón  y  aquellas  que  involucran  especies 
altamente  reactivas,  frecuentemente  dependen  de  la  difusión  [31]  y  ocurren  en  procesos 
con  barreras  de  activación  pequeñas  que  hacen  que  la  velocidad  se  limite  por  la  rapidez 
con  la  que  se  difunden  los  reactivos  a  través  del  disolvente  [50];  por  ello,  se  empleó  la 
corrección  por  el  límite  de  difusión  de  Collins-Kimball  [51]  según  la  Ec.  (8), 

.  _  kTST  kD 

KaV  ~  kTST+kD  (8> 

donde  kap  es  la  constante  de  velocidad  aparente;  l^ST  es  la  constante  de  velocidad  según 
la  TST  y  kD  es  la  constante  de  velocidad  de  Smoluchowski  o  constante  de  difusión;  esta 
se  calcula  por  medio  de  la  Ec.  (9), 

kD  -  4 kRabDabNa  (9) 

donde  RAB  es  la  distancia  de  reacción  que  se  obtiene  a  partir  del  TS,  al  medir  la  distancia 
entre  los  átomos  que  transfieren  el  H  y  la  distancia  entre  aquellos  que  forman  el  aducto 
radical  para  las  reacciones  que  siguen  el  mecanismo  HT  y  RAF,  respectivamente.  Para 

12 

Santiago  Alvarez  (2019) 


el  mecanismo  SET,  la  distancia  de  reacción  se  obtiene  a  partir  de  la  suma  de  los  radios 
de  los  reactivos  participantes.  NA  es  el  número  de  Avogadro  y  DAB  es  el  coeficiente  de 
difusión  mutua  de  los  reactivos  A  y  B.  DAB=DA+DB  se  calcula  por  medio  de  la  Ec.  (10)  de 
Stokes-Einsteln, 


D 


AoB 


kJ 

6^WaoB 


(10) 


donde  r|  es  la  viscosidad  del  disolvente,  en  este  caso  t|  =  8,91  x  1 0-4  Pa.s  para  el  agua  y 
rj  a*  8,62  x  10~4  Pa.s  para  el  etanoato  de  pentilo;  aAoB  representa  el  radio  del  soluto  y  las 
unidades  de  DAoB  y  kD  son  m2  s_1  y  L  mol-1  s'1,  respectivamente. 

En  Eyringpy  se  encuentran  implementadas  las  Ec.  (5)  -  (10)  y  las  demás  funciones  de 
la  Termodinámica  Estadística,  necesarias  para  calcular  las  constantes  de  velocidad. 
En  el  archivo  de  entrada  de  Eyringpy  se  ajustan  el  método,  el  efecto  túnel,  la  difusión, 
temperatura,  pH,  distancias  de  reacción  y  disolvente,  junto  con  los  archivos  de  salida  (out) 
correspondientes  a  los  reactivos,  productos  y  estado  de  transición  que  arroja  Goussian. 
Si  todos  los  valores  de  entrada  son  correctos  y  sus  cálculos  convergen,  entonces  el 
programa  arroja  un  archivo  de  salida  con  los  correspondientes  datos  termoquímicos 
y  cinéticos. 


RESULTADOS 

El  3,4-dlhidroxibenzaldehído  se  representa  con  la  notación  EI2C7H403  y  sus  especies 
ionizadas  correspondientes  como  EHC7H403_  y  C7H4032~.  El  resultado  de  la  búsqueda 
conformaclonal  indica  que  las  moléculas  que  forman  puentes  de  hidrógeno 
intramoleculares  y  tienen  el  oxígeno  del  grupo  aldehido  en  la  misma  orientación  que 
los  OH  en  el  plano  del  anillo  aromático,  son  estructuras  de  mínima  energía. 

Para  el  H2C7H403,  el  oxígeno  de  la  posición  3  forma  un  puente  de  hidrógeno  de  2,155 
Á  que  estabiliza  al  sistema.  Cuando  el  H2C7H403  se  disocia,  el  hidrógeno  del  OH  de 
la  posición  4  se  desprende  para  formar  el  HC7H403",  la  distancia  C-0  se  contrae  0,08 
Á  respecto  del  3,4-dihidroxibenzaldehído  sin  disociar  y  el  hidrógeno  del  OH  de  la 
posición  3  gira  hacia  el  oxígeno  adyacente  para  formar  un  puente  de  hidrógeno  de 
2,045  Á.  En  la  Fig.  3  se  observa  que  cuando  el  3,4-dlhidroxibenzaldehído  pierde  los 
hidrógenos,  las  distancias  de  enlace  C-0  de  las  posiciones  3  y  4  disminuyen  0,05  Á  y  0,07 
Á  respectivamente.  Para  el  C7H4032_  el  isómero  de  mínima  energía  es  aquel  cuyo  átomo 
de  oxígeno  del  grupo  aldehido  tiene  la  misma  orientación  del  oxígeno  en  la  posición 
3.  Cada  estructura  de  la  figura  hace  las  veces  de  reactivo  en  los  diferentes  mecanismos 
que  se  estudian  en  esta  investigación.  Los  cambios  y  la  reactividad  de  estas  moléculas 
se  indican  en  la  Fig.  3. 
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FIGURA  3.  Estructura  de  las  especies  del  3,4-dihidroxibenzaldehído  con  el  nivel  M05-2X/6-31 1  +G(d,p).  Las 
distancias  de  enlace  se  expresan  en  Á.  El  1  denota  al  H2C7H403,  el  2  al  HC7H403"  y  el  3  al  C7H4032'.  Los  valores  entre 
paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 


Mecanismo  de  transferencia  de  hidrógeno 

El  3,4-dihidroxibenzaldehído  tiene  seis  hidrógenos  en  su  estructura,  los  cuales  pueden 
ser  removidos  por  el  HOO.  Las  estructuras  de  losTS  para  este  mecanismo  se  muestran  en 
la  Fig.  4.  En  etanoato  de  pentilo,  los  átomos  de  hidrógeno  del  3,4-dihidroxibenzaldehído 
a  transferirse  se  ubican  más  cerca  del  átomo  de  oxígeno  del  radical  hidroperoxilo. 
Para  extraer  el  hidrógeno  del  grupo  aldehido  y  los  hidrógenos  2,  5  y  6  enlazados  a 
los  carbonos,  el  HOO-  se  coloca  junto  al  plano  del  anillo  aromático;  mientras  que  para 
extraer  los  hidrógenos  de  los  grupos  OH  de  las  posiciones  3  y  4,  el  radical  se  dispone 
sobre  el  plano  del  anillo. 

Los  estados  de  transición  correspondientes  a  las  reacciones  de  transferencia  del 
hidrógeno  del  grupo  aldehido  y  los  hidrógenos  de  los  OH  de  las  posiciones  3  y  4  son 
los  de  energía  más  baja.  Los  hidrógenos  de  los  grupos  OH  son  más  ácidos  que  el  H  del 
grupo  aldehido  y,  por  tanto,  su  transferencia  es  más  fácil.  Ahora,  entre  el  H3  y  el  H4, 
cuando  el  HOO-  se  aproxima,  el  H4  gira  y  se  rompe  el  puente  de  hidrógeno  formado  con 
el  oxígeno  del  OH  de  la  posición  3  y,  por  ello,  la  transferencia  del  H4  cuesta  2,6  kcal/mol 
más  que  la  del  H3  en  medio  acuoso  y  4,0  kcal/mol  más  en  medio  no  polar. 

Los  hidrógenos  2, 5  y  6  enlazados  a  los  átomos  de  carbono  no  son  ácidos;  su  transferencia 
presenta  barreras  de  energía  mayor  a  las  20  kcal/mol  y,  por  ello  ,  un  medio  de  reacción 
fisiológico  con  temperatura  aproximada  de  310  K  no  proveería  tal  demanda  energética. 
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TS1-H3  6.3  (3.6) 


TS1-H  8.9  (7.2) 


TS1-H4  8.9  (7.6) 


TS1-H6  25.2(23.2) 


TS1-H5  24.6(24.8) 


TS1-H2  25.2(25.9) 


FIGURA  4.  Estructura  de  los  TS  que  se  forman  entre  1  y  el  H00-  en  el  mecanismo  HT  con  el  nivel  M05-2X/6- 
31 1  +G(d,p).  Las  distancias  en  Á.  En  la  parte  inferior  de  cada  estructura  se  indica  el  valor  de  la  energía  en  kcal/mol 
con  respecto  a  los  reactivos.  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 


El  parámetro  termodinámico  que  determina  la  viabilidad  de  los  mecanismos  HT  es  el 
cambio  de  energía  libre  de  Gibbs;  si  esta  es  menor  a  cero,  la  reacción  es  exergónica;  si  es 
cero,  la  reacción  se  encuentra  en  equilibrio,  y,  si  es  positiva,  la  reacción  es  endergónica.  En 
la  Tabla  1  se  muestran  los  valores  termoquímicos  para  los  procesos  HT  entre  1  y  el  HOO- 
en  ambos  disolventes.  La  energía  libre  de  Gibbs  de  los  complejos  de  reactivos  (AGCR)  es 
menor  que  la  de  los  reactivos.  Esto  ocurre  porque  las  moléculas  forman  interacciones 
intramoleculares  fuertes  entre  el  aldehido  y  el  radical  que  varían  entre  7  y  9  kcal/mol  y 
estabilizan  el  complejo.  En  cambio,  la  energía  de  Gibbs  de  los  complejos  de  productos 
(AGCP)  es  en  general  mayor,  debido  a  que  el  aumento  de  la  entropía  en  el  rompimiento 
de  enlaces  en  los  procesos  HT,  no  compensa  la  disminución  de  la  energía  electrónica 
corregida  con  la  energía  de  punto  cero,  de  acuerdo  con  la  relación  AG=AH  -  TAS. 
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TABLA  1.  Cambio  de  energía  de  Gibbs  en  kcal  mol'1  para  el  mecanismo  HT  ente  1  y  el  H00-  a  298,15  K  * 


AG 

CR 

AG* 

TS 

AG 

CP 

AGp  1 

H 

-9,7  (-13,6) 

8,9  (7,2) 

-0,8  (-3,1) 

-7,3  (-9,3) 

H2 

-6,0  (-9,2) 

25,2  (25,7) 

20,8(18,9) 

15,4(15,7) 

H3 

-8,3  (-9,2) 

6,3  (3,6) 

-9,9  (-12,7) 

-15,0  (-17,0) 

H4 

-6,9  (-10,6) 

8,9  (7,6) 

-7,5  (-7,5) 

-12,8  (-12,0) 

H5 

-8,3  (-10,0) 

24,6  (24,8) 

22,0  (22,1) 

15,0(15,5) 

H6 

-7,1  (-10,2) 

23,7  (23,2) 

20,5  (21,0) 

14,1  (13,9) 

*  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 


La  energía  de  Gibbs  de  reacción,  calculada  como  la  diferencia  entre  los  reactivos  y 
productos  (AGP)  indica  que  se  favorecen  las  reacciones  HT  en  el  H  del  grupo  aldehido  y 
la  de  los  hidrógenos  de  los  OH  de  las  posiciones  3  y  4  en  ambos  disolventes.  El  proceso 
más  exergónico  y  que  presenta  la  energía  de  activación  más  baja,  es  la  reacción  de 
transferencia  del  H3.  La  energía  de  activación  efectiva  (calculada  desde  el  complejo 
de  reactivos)  es  14,6  kcal/mol  para  la  reacción  en  medio  acuoso.  En  la  tabla,  H  hace 
referencia  al  hidrógeno  del  grupo  aldehido;  H2,  H5,  H6  a  los  hidrógenos  del  anillo 
aromático;  H3  y  H4  a  los  hidrógenos  de  los  grupos  OH  en  las  posiciones  3  y  4. 


En  la  Fig.  5  se  muestra  el  perfil  de  reacción  de  transferencia  del  H3  en  medio  acuoso.  El 
complejo  de  reactivos  se  forma  por  medio  de  un  puente  de  hidrógeno  entre  el  H  del 
radical  hidroperoxilo  y  el  O  en  el  grupo  aldehido;  como  productos  se  obtiene  el  radical 
en  el  oxígeno  de  la  posición  3  más  peróxido  de  hidrógeno.  Termodinámicamente,  el 
complejo  se  encuentra  a  -8,3  kcal/mol  con  respecto  a  los  reactivos,  entonces  la  energía 
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de  activación  efectiva  se  calcula  desde  el  complejo  hasta  el  TS  que  se  encuentra  a  6,3 
kcal/mol.  La  barrera  energética  real  de  este  proceso  es  14,6  kcal/mol:  la  más  baja  de 
todos  los  procesos  HT,  y  con  la  reacción  más  negativa  en  la  energía  de  Gibbs.  Se  espera 
que  las  reacciones  que  presentaron  energía  de  Gibbs  más  negativa,  originen  productos 
mayoritarios.  El  tamaño  de  las  barreras  de  energía  y  la  formación  de  complejos  influye  en 
la  rapidez  de  las  reacciones,  como  se  indica  a  continuación. 

Las  constantes  de  velocidad  para  las  reacciones  de  transferencia  de  hidrógeno  entre  1 
y  el  HOO-  se  muestran  en  la  Tabla  2.  El  parámetro  indica  la  constante  de  velocidad, 
sin  considerar  ninguna  corrección.  En  los  casos  donde  existe  tunelamiento,  el  factor  k  se 
multiplica  por  esta  constante  que  da  lugar  a  .  Los  valores  de  kap  consideran  la  corrección 
por  difusión  y  tunelamiento,  en  caso  de  existirlo. 

TABLA  2.  Constantes  de  velocidad  para  el  mecanismo  de  transferencia  de  hidrógeno  en  IVL1  s_1  para  1  frente  al 

HOO- a  298,15  K.* 


fcTST 

K 

eck 

ktst 

K 

k 

ap  J 

H 

4,6x1 07  (8,4x1 08) 

24,0  (20,6) 

1,1  x  109(1,7x  10'°) 

2,8x1 09  (2,9  x1o9) 

7,9x  1 08  (2,5x1 09) 

H2 

4,9x1 0"5  (1,7x1  a5) 

2,8x1 09  (2,9  x1o9) 

4,9x10"  (1,7x10") 

H3 

3,7x1 09  (3,3x10") 

413,1  (114,9) 

1,5  x  10"  (3,8x1 013) 

2,6x1 09  (2,7  x1o9) 

2,6x  1 09  (2,7x1 09) 

H4 

4,8x1 07  (3,9x  1 08) 

3,0x  1 03  (165,1) 

1,5x10"  (6,4x1 010) 

2,7x1 09  (2,7  x1o9) 

2,6x  1 09  (2,6x1 09) 

H5 

1,3x10"  (1,1x10") 

2,8x1 09  (2,9  x1o9) 

1,3  x  10"(1,1  x  10") 

H6 

6,2x10"  (1,5x1 0'3) 

2,8x1 09  (2,9  x1o9) 

6,2  x  10"  (1,5  x  10") 

*  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 


En  los  resultados  termoquímicos  se  muestra  que  los  procesos  HT  de  los  hidrógenos  ácidos 
en  medio  de  etanoato  de  pentilo  presentan  valores  de  energía  de  activación  más  baja  que 
en  medio  acuoso;  esto  hace  que  el  valor  de  sea  más  alto  en  medio  no  polar.  En  ambos 
medios,  la  reacción  en  H3  es  tan  rápida  que  su  constante  tiende  al  límite  de  difusión  con 
una  rapidez  de  2,6  x  109  IVL1  s_1  para  agua  y  2,7  x  109  M_1  s_1  para  etanoato  de  pentilo.  En  la 
reacción  de  transferencia  del  hidrógeno  del  OH  de  la  posición  4,  el  efecto  túnel  que  se 
origina  por  la  formación  de  los  complejos  de  reactivos  da  lugar  a  valores  altos  de  k;  aumenta 
la  rapidez  de  reacción  y  produce  que  la  constante  de  velocidad  tienda  a  la  de  difusión. 

Las  reacciones  de  transferencia  de  los  hidrógenos  2,  5  y  6  presentan  energías  de 
activación  elevadas;  por  ello,  las  constantes  de  velocidad  tienden  a  cero  y  no  tienen 
relevancia  en  la  reactividad  de  ninguna  de  las  especies  del  3,4-dihidroxibenzaldehído 
frente  al  HOO*. 

De  Igual  forma,  para  2  y  el  radical  hidroperoxilo,  los  caminos  de  reacción  más  favorecidos 
son  la  transferencia  del  hidrógeno  del  grupo  aldehido  y  el  hidrógeno  del  OH  de  la  posición 
3  en  ambos  disolventes.  La  energía  de  activación  efectiva  más  baja  es  de  10,7  kcal/mol 
en  medio  acuoso  (Fig.  6)  mientras  que  en  medio  apolar  es  de  1 8,9  kcal/mol:  ambas  para 
el  proceso  HT  en  el  H3.  La  barrera  para  H3,  en  medio  no  polar,  es  significativamente  más 
grande  que  en  medio  polar;  esto  explica  por  qué  es  más  difícil  que  ocurra  la  reactividad 
de  especies  ionizadas  en  disolventes  apolares. 
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FIGURA  6.  Perfil  de  reacción  para  2  siguiendo  el  mecanismo  HT  en  medio  acuoso. 

Para  el  caso  de  3,  la  reacción  viable  es  la  transferencia  del  H  en  el  grupo  aldehido,  la 
energía  de  activación  efectiva  es  de  1 7,7  y  1 9,6  kcal/mol  con  constantes  de  velocidad  de 
2,0  x  1 09  y  2,8  x  1 09  M_1  s_1  en  medio  acuoso  y  no  polar,  respectivamente. 

Mecanismo  de  formación  de  aductos 

El  3,4-dihidroxibenzaldehído  posee  un  átomo  de  carbono  en  el  grupo  aldehido  y  6  en  el 
anillo  aromático  en  donde  puede  darse  la  adición  del  radical  hidroperoxilo.  La  estructura 
de  losTS  se  muestra  en  la  Fig.  7.  En  las  reacciones  de  formación  de  los  aductos  de  los 
carbonos  1 , 3  y  4,  el  radical  HOO-  se  adiciona  de  tal  manera  que  el  grupo  aldehido  y  los 
grupos  hidroxilo  de  esas  posiciones  se  inclinan  un  ángulo  aproximado  de  15°,  debido  a 
las  repulsiones  electrónicas  entre  los  oxígenos. 

La  reacción  de  formación  del  aducto  en  la  posición  4  presenta  el  estado  de  transición 
TS1-C4  como  el  de  menor  energía,  le  sigue  el  TS1-C6  correspondiente  a  la  adición 
del  radical  en  el  C6.  Esta  tendencia  indica  que  el  mecanismo  RAF  se  favorece  en  las 
posiciones  orto  y  poro  en  el  anillo  aromático,  respecto  del  grupo  aldehido.  Según  este 
análisis, el  siguiente  estado  de  transición  sería  elTS1-C2en  lugar  del  TS1 -C3;  sin  embargo, 
la  adición  del  radical  en  el  C3  hace  que  el  anillo  aromático  no  quede  totalmente  plano, 
por  lo  que  la  distancia  entre  el  oxígeno  de  la  posición  4  y  el  hidrógeno  del  OH  de  la 
posición  3  aumenta  a  2,1 95  Á,  al  romper  el  puente  de  hidrógeno  y,  por  ello,  existe  una 
diferencia  de  0,7  kcal/mol  entre  el  TS1 -C3  y  el  TS1 -C2.  El  TS,  correspondiente  a  la  adición 
en  el  C5  y  los  demás  procesos  RAF,  presenta  energía  más  elevada  y  dificulta  la  formación 
de  los  productos  en  condiciones  fisiológicas. 
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TS1-C4  9.6(10.4) 


TS1-C3  12.5(12.6) 


TS1-C5  14.6(15.7) 


TS1-C6  11.3(11.9) 


TS1-C2  13.2(13.4) 


TS1-C1  15.2(14.8) 


TS1-C  20.4(20.5) 

FIGURA  7.  Estructura  de  los  TS  que  se  forman  entre  1  y  el  HOO-  en  el  mecanismo  RAF  con  el  nivel  M05-2X/6- 
31 1+G(d,p).  Las  distancias  en  Á.  En  la  parte  inferior  de  cada  estructura  se  indica  el  valor  de  la  energía  en  kcal/mol 
con  respecto  a  los  reactivos.  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 
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De  acuerdo  a  la  Tabla  3,  los  complejos  de  reactivos  se  favorecen  en  todos  los  procesos  de 
1  adición  y  los  de  productos  -en  este  caso-  no  se  forman,  ya  que  se  trata  de  productos  de 

encfenáfe  adición  y  no  de  moléculas  separadas  como  en  el  caso  del  mecanismo  HT,  donde  había 

m9emenas  posibilidad  de  formar  puentes  intermoleculares.  El  cambio  de  energía  de  Gibbs  entre 
productos  y  reactivos  indica  que  la  formación  del  aducto  en  el  carbono  4  (posición  paro) 
es  exergónica  con  un  AGp  =  -0,4  kcal/mol  en  medio  acuoso.  En  comparación  con  los 
productos  del  mecanismo  HT,  los  de  adición  no  son  mayoritarios,  ya  que  la  energía  de  Gibbs 
es  ligeramente  negativa  y,  por  tanto,  la  constante  de  equilibrio  (K)  es  pequeña  según  la 
relación  AG  =  -  RT  Ln(K).  En  la  tabla,  el  literal  C  hace  referencia  al  carbono  del  grupo  aldehido. 


TABLA  3.  Cambio  de  energía  de  Gibbs  en  kcal  mol'1  para  el  mecanismo  RAF  entre  1  y  el  H00-  a  298,15  K.* 


AG 

CR 

AG* 

TS 

AGp 

c 

-10,0  (-14,5) 

20,4(20,5) 

14,7  (8,1) 

C1 

-6,6  (-9,7) 

15,2(14,8) 

10,8  (9,7) 

C2 

-7,0  (-13,3) 

13,2(13,4) 

4,6  (0,2) 

C3 

-7,4  (-10,2) 

12,5(12,6) 

4,0  (3,1) 

C4 

-7,1  (-9,6) 

9,6  (10,4) 

-0,4  (0,4) 

C5 

-6,2  (-10,4) 

14,6(15,7) 

8,5  (8,5) 

C6 

-6,6  (-10,4) 

(11,9) 

2,2  (3,0) 

*  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 

La  Tabla  4  indica  las  constantes  de  velocidad  para  el  mecanismo  de  formación  de  aductos 
para  cada  carbono  de  1,  frente  al  radical  hidroperoxilo.  En  ninguno  de  estos  procesos 
existe  tunelamiento:  el  HOO-  no  puede  atravesar  las  barreras  de  energía  debido  a  la 
masa  relativamente  alta  que  posee.  En  general,  las  reacciones  de  formación  de  aductos 
no  son  espontáneas  y,  si  ocurriesen,  sus  constantes  de  velocidad  son  bajas  o  próximas  a 
la  tasa  de  extracción  de  H  en  los  ácidos  grasos  presentes  en  las  células. 

TABLA  4.  Constantes  de  velocidad  (M1  s1)  para  el  mecanismo  RAF  de  1  frente  al  H00-  a  298,15  K* 


kTST 

ko 

* 

ap 

C 

1,7x1o'1  (1,3  x1o'1) 

2,0  x  1 09  (1,9  x  109) 

17x1 0-1  (1,3x1 0-1) 

C1 

1,1  xIO3  (2,3  xIO3) 

2,0 x  1 09  (2,1  xIO9) 

1,1  x  1 03  (2,3  x  1 03) 

C2 

3,1  xl  04  (2,5x  1 04) 

2,1  x  109  (2,1  x  109) 

3,1  x  1 04  (2,5  x  1 04) 

C3 

1 ,0 x  1 05  (8,3 xIO4) 

2,1  xIO9  (2,2 x  1 09) 

1 ,0  x  1 05  (8,3  x  1 04) 

C4 

1 ,3  x  1 07  (3,4  x  1 06) 

2,1  x  109  (2,2  x  1 09) 

1 ,3  x  1 07  (3,4  x  1 06) 

C5 

2,9x1 03  (4,5x1 02) 

2,1  x  109  (2,1  xIO9) 

2,9x1 03  (4,5x1 02) 

C6 

7,9x  1 05  (3,2x  1 05) 

2,1  xIO9  (2,2 x  1 09) 

7,9x  1 05  (3,2x  1 05) 

*  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 


Los  procesos  RAF  para  2  presentan  mayor  número  de  reacciones  con  energía  de  Gibbs 
negativa:  esto  se  debe  a  que  el  puente  de  hidrógeno  que  estabiliza  las  estructuras  se 
mantiene  en  los  productos.  En  este  caso,  la  adición  en  el  carbono  del  grupo  aldehido, 
los  carbonos  2  y  3  en  medio  acuoso,  al  igual  que  los  carbonos  2,  3,  4  y  6  en  medio  no 
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polar,  son  exergónicos.  Cuando  el  disolvente  es  etanoato  de  pentllo,  los  complejos  de 
reactivos  presentan  energía  más  baja  en  comparación  con  los  del  medio  acuoso.  En 
la  Fig.  8  se  muestra  la  reacción  de  formación  del  aducto  radical  en  el  carbono  2;  esta 
es  exergónica  para  ambos  disolventes,  con  un  cambio  de  energía  de  Gibbs  de  -0,4  en 
agua  y  -3,7  kcal/mol  en  medio  apolar.  La  energía  de  activación  efectiva  es  16,2  kcal/ 
mol  en  medio  acuoso  y  21,3  kcal/mol  en  medio  no  polar;  debido  a  que  el  complejo  de 
reactivos  en  medio  apolar  se  encuentra  a  aproximadamente  5,4  kcal/mol  más  abajo  en 
la  superficie  de  energía  potencial.  Estos  datos  termoquímicos  explican,  además,  por  qué 
las  reacciones  de  especies  ionizadas  en  disolventes  apolares  no  ocurren. 


Las  reacciones  de  adición  en  el  C7H4032_  presentan  barreras  de  energía  que  tienden  a  0 
kcal/mol  con  respecto  a  3  +  HOO-. 


Complejo  2(C2) 
-8.2  (-13.6) 


TS2-C2  2A(C2) 

8.0  (7.7)  -0.4  (-3.7) 


FIGURA  8.  Estructura  del  complejo  de  reactivos,  estado  de  transición  y  producto  de  adición  correspondientes  a  la 
reacción  deformación  del  aducto  en  C2  con  el  nivel  M05-2X/6-311+G(d,p).  En  la  parte  inferior  de  cada  estructura 
se  indica  el  valor  de  la  energía  en  kcal/mol  con  respecto  a  2  +  H00-.  Los  valores  entre  paréntesis  corresponden  al 

medio  de  etanoato  de  pentilo. 

Mecanismo  de  transferencia  de  electrones 

El  H2C7EI403  no  participa  en  el  mecanismo  de  transferencia  de  electrones;  estos 
procesos  resultan  demasiado  endergónicos,  ya  que  es  difícil  transferir  electrones  desde 
una  especie  neutra  como  el  H2C7H403  hacia  el  radical.  El  valor  de  energía  de  Gibbs  de 
reacción  es  23,8  kcal/mol  cuando  el  disolvente  es  agua  y  67,4  kcal/mol  en  medio  de 
etanoato  de  pentilo.  La  energía  de  activación  que  deben  alcanzar  estos  procesos  es  23,9 
y  81 ,6  kcal/mol,  en  medio  acuoso  y  no  polar,  respectivamente. 


La  termoquímica  de  los  procesos  SET  refleja  la  dificultad  que  presentan  las  especies 
ionizadas  para  reaccionar  en  disolventes  apolares;  así,  para  2  en  medio  acuoso  la 
reacción  (1 1)  es  exergónica  con  un  AG=  -2,5  kcal/mol  y  una  energía  de  activación  de 
4,5  kcal/mol;  mientras  que  en  medio  no  polar,  la  reacción  es  endergónica  con  AG=  1 1,8 
kcal/mol  y  una  barrera  de  13,4  kcal/mol. 
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avances 
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En  general,  3  presenta  valores  elevados  de  las  constantes  de  transferencia  de 
electrón  que  en  principio  hacen  pensar  que  el  C7H4032_  es  la  especie  predominante 
en  el  sistema;  sin  embargo,  al  analizar  el  diagrama  de  distribución  de  especies  del 
3,4-dihidroxibenzaldehído  (Fig.9)  se  observa  que  en  condiciones  fisiológicas  (pH=7,4)  el 
C7H4032_  tiene  una  fracción  molar  de  1 ,8  x  1 0-5'  esto  significa  que  la  concentración  de  esta 
especie  y  las  reacciones  que  genera  son  prácticamente  despreciables. 


FIGURA  9.  Diagrama  de  distribución  de  especies  del  3,4-dihidroxibenzaldehído  en  función  del  pH. 

Actividad  antioxidante  del  3,4-dihidroxibenzaldehído 

La  fracción  molar  del  HC7H403_  es  0,602  y  la  del  H2C7H403  es  0,398,  lo  que  indica  que  el 
monoanión  es  la  especie  que  se  encuentra  en  mayor  proporción.  Ahora,  al  considerar 
todos  los  mecanismos  se  deduce  que  el  mecanismo  HT  es  el  principal  camino  de 
reacción  que  tomaría  el  3,4-dihidroxibenzaldehído  para  reducir  la  toxicidad  del  radical 
hidroperoxilo,  ya  que  estos  procesos  reportan  las  constantes  de  velocidad  con  los 
órdenes  de  magnitud  más  altos,  los  valores  de  energía  de  activación  efectiva  más  bajos, 
así  como  los  valores  de  cambio  de  energía  libre  de  Gibbs  más  negativos. 

Para  el  proceso  en  medio  acuoso,  en  el  cual  el  monoanión  transfiere  un  protón  al  radical 
hidroperoxilo,  se  tiene  una  constante  de  velocidad  k1  según  la  reacción  (1 2). 

HC.Hfi;  +H00»^^C1H,0^+H202  (i2) 

En  una  reacción  paralela,  el  HC7H403"  transfiere  un  electrón  hacia  el  HOO-  con  una 
constante  de  velocidad  k2  según  la  reacción  (1 3). 

HC.Hfi;  +  HOO»^^HC1HaOj»+HOO~  (i  b) 

Adicionalmente,  la  especie  C7H403-~,  producto  de  la  transferencia  de  hidrógeno  también 
cede  un  electrón  al  radical  según  la  reacción  (14). 

C7//403  •“  +HOO* — » C7tf403  +HOO~  (14) 
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Las  ecuaciones  diferenciales  que  gobiernan  estos  procesos  se  reflejan  en  las  Ec.  (1 5)  y  (1 6) 

d[H°°‘]  =  — Ar, [HC7H403~  ] [HOO‘] - k2 \HCnHtO{ ] [HOO*\ -k^Cpp,-  ][HOO-]  (1  5) 
dt 

d\£jH£¿_  ]  =  ki[HC7HA0;][H00-]-k3[C1H403'-][H00.].  (16) 
dt 

Para  encontrar  [C7H403-_]  se  emplea  la  aproximación  del  estado  estacionario  en  la  Ec.  (1 5): 

k^HC.Hfi^HOO*}  =  k^HfipWiOO'}  (i  7) 


entonces: 


[c7//4o3.-]= 


k^HC.H.o;] 


(18) 


al  reemplazar  la  expresión  de  la  Ec.  (1 8)  en  la  Ec.  diferencial  (1 6)  se  obtiene  la  Ec.  (1 9): 


d^HO°^  =  +k1)[HC1H40;][H00>]-k,  ^HCiH4°i  }[HOO>]  (19) 

dt  k3 

y  al  simplificar  se  obtiene  la  Ec.  (20): 

=  -(2kt  +  k2)[HC7H4Oi][HOOp  (20) 

dt 

Como  puede  notarse  en  la  Ec.  (20),  la  constante  de  velocidad  global  para  la  inhibición 
del  radical  HOO-  corresponde  a  +  k2.  El  valor  de  k1  es  2,8  x  109  M_1  s_1  mientras  que 
el  valor  de  k2es  2,2  x  109  M_1  s_1  (ver  Tabla  S.2  en  Material  Suplementario);  es  decir,  la 
constante  global  tiene  un  valor  de  7,8  x  109  M_1  sT  Al  constituir  la  Ec.  (20)  se  obtiene  la 
ley  integrada  de  velocidad  del  proceso  de  acuerdo  con  la  Ec.  (21). 

=  (2 K  +k2)([HOO-\  ~[HC7H4O3-]0)t.  (21) 


donde  [HOO-]0  y  [HC7H4O3_]0son  las  concentraciones  iniciales  del  radical  hidroperoxilo 
y  el  monoanión,  respectivamente,  siendo  [HC7H4O3_]0=  0,602CH2C7H4O3  (ver  Material 
Suplementario).  Esta  ecuación  permite  conocer  las  concentraciones  instantáneas  del 
radical  y  el  aldehido  en  fundón  del  tiempo. 

En  ambiente  no  polar,  el  H2C7H403es  la  especie  química  que  predomina  al  transferir  un 
hidrógeno  al  radical  según  la  reacción  (22). 

H2C7H403  +  HOO-^  HC7H403  .+tf202  (22) 

Esta  reacción  tiene  una  constante  de  velocidad  k  =  2,7  x  109M_1  s_1  (Tabla  2);  al  considerar 
que  las  constantes  de  velocidad  correspondientes  al  daño  que  causa  el  HOO-  en  los 
ácidos  grasos  poliinsaturados,  están  en  el  rango  de  1  -3  x  1 03  M_1  s_1  [1 2],  se  predice  que 
el  3,4-dlhidroxibenzaldehído  tiene  la  capacidad  de  inhibir  estos  radicales. 


Ln 


[//C7//46>3-]0[MX>] 

\HOO*UHC,H.On 
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3CÍ  DISCUSIÓN 

encfen cíale  La  capacidad  antioxidante  del  3,4-dihidroxibenzaldehído  frente  al  radical  hidroperoxilo  se 

ingenierías  estuc|¡Q  p0r  medio  de  un  análisis  termoquímico  y  cinético  en  dos  disolventes:  agua  para 
imitar  el  ambiente  acuoso  del  citosol  presente  dentro  de  la  células  y  etanoato  de  pentilo 
para  imitar  el  medio  no  polar  de  la  membrana  celular,  donde  se  generan  estos  radicales. 

El  análisis  conformacional  indica  que  las  especies  químicas  más  estables  para  el  3,4-dihi¬ 
droxibenzaldehído  poseen  estructuras  planas  con  un  enlace  de  hidrógeno  intramolecular. 
Este  resultado  coincide  con  un  análisis  que  se  realizó  con  el  nivel  de  teoría  B3LYP/6-31G  [52]. 

Termodinámicamente,  se  encontró  que  las  reacciones  de  transferencia  de  los 
hidrógenos  ácidos  del  compuesto  fenólico  presentan  los  cambios  de  energía  libre  de 
Gibbs  de  reacción  más  negativos  y  las  barreras  de  energía  más  bajas,  lo  que  significa 
que  el  mecanismo  HT  es  la  vía  principal  que  toma  la  reacción.  A  medida  que  aumenta 
el  grado  de  deprotonación  del  aldehido,  se  favorecen  las  reacciones  de  adición  en  los 
átomos  de  carbono  de  la  molécula;  sin  embargo,  la  exergonicidad  de  estos  procesos 
sigue  siendo  menor  a  los  procesos  HT. 

La  constante  de  velocidad  para  la  reacción  entre  el  ácido  3,4-dihidroxibenzóico  y  el 
HOO-  es  1 ,26  x  1 07  M_1  s_1  en  medio  acuoso  y  5,1 4  x  1 03  M_1  s_1  en  medio  de  etanoato  de 
pentilo  [39]:  esto  sugiere  que  la  capacidad  antioxidante  del  3,4-dihidroxibenzaldehído  es 
mayor  que  la  del  ácido  3,4-dihidroxibenzóico.  Los  ensayos  experimentales  de  capacidad 
antioxidante  equivalente  de  trolox  (TEAC)  y  el  ensayo  de  2,2-difenM-l-picrililhidracilo 
(DPPH-)  [53]  junto  con  los  cálculos  de  varios  descriptores  moleculares  en  DFT  realizados 
para  estos  compuestos,  confirman  esta  tendencia  [52].  Los  resultados  en  ambiente 
no  polar  concuerdan  con  las  propiedades  de  los  antioxidantes  fenólicos  de  retardar  la 
oxidación  de  lípidos  a  través  del  mecanismo  de  transferencia  de  hidrógeno  [54]. 

El  mecanismo  secuencial  de  transferencia  de  los  hidrógenos  3  y  4  del  aldehido  se 
estudia,  además,  desde  el  punto  de  vista  teórico;  siempre  y  cuando  se  considere  el 
carácter  multirreferencial  del  proceso:  su  importancia  radica  en  la  determinación  de  la 
capacidad  prooxidante  del  compuesto. 
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MATERIAL  SUPLEMENTARIO 


Tabla  S.1 .  Cambio  de  energía  de  Gibbs  en  kcal  mol-1  para  2  frente  al  H00-  a  298,15  K .* 


Mecanismo 

AG 

CR 

AG* 

TS 

AG 

CP 

Q_ 

< 

HT 

H 

-11,8  (-19,5) 

7,2  (2,7) 

-2,5  (-6,9) 

-8,1  (-11,3) 

H2 

-11,3  (-15,6) 

20,6(18,8) 

16,8(12,3) 

13,1  (12,5) 

H3 

-8,0  (-16,2) 

2,7  (2,7) 

-1 7,2  (-24,8) 

-21,3  (-19,8) 

H5 

-11,4  (-13,5) 

20,4(17,9) 

17,8(10,1) 

13,1  (12,4) 

H6 

-11,2  (-14,0) 

21,1  (17,3) 

18,5(13,6) 

11,3(9,9) 

RAF 

C 

-8,5  (-13,5) 

4,2  (6,1) 

-4,7  (0,9) 

C1 

-7,7  (-13,6) 

10,1  (8,0) 

8,6  (6,5) 

C2 

-8,2  (-13,6) 

8,0  (7,7) 

-0,4  (-3,7) 

C3 

-8,4  (-13,4) 

4,0  (4,0) 

-1,1  (-3,1) 

C4 

-8,2  (-12,9) 

6,0  (8,0) 

3,8  (-7,9) 

C5 

-7,7  (-1 3,4) 

7,7  (6,5) 

4,4  (2,4) 

C6 

-8,3  (-13,5) 

9,4  (7,3) 

2,1  (-0,7) 

SET 

4,5(13,4) 

-2,5(11,8) 

*  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 


Tabla  S.2.  Constantes  de  velocidad  (M1  s-1)  para  2  frente  al  H00-  a  298,15  K  * 


kTST 

K 

1,  eck 

ktst 

K 

k 

ap 

HT 

H 

7,4x108  (1,7x1 012) 

12,6  (6,1) 

9,4x1 09  (1,0x1 013) 

2,9x1 09 

(2,8  x1o9) 

2,2x  1 09  (2,8x  1 09) 

H2 

1,2  x1o-1  (2,7) 

09,5) 

(5,2  x1o1) 

2,9x1 09 

(2,8x1 09) 

1,2  x1o-1  (5,2  x1o1) 

H3 

1 ,5  x  1 012  (1 ,6  x  1012) 

1 ,2  x  1 04  (11,6) 

1 ,9  x  1 016  (1 ,9  x  1 013) 

2,8x1 09 

(2,8  x1o9) 

2,8x  1 09  (2,8x  1 09) 

H5 

1 ,6  x  1 01  (12) 

(17,7) 

(2,2x1 02) 

2,9x1 09 

(2,8  x1o9) 

1,6  x1o-1  (2,2x1 02) 

H6 

5,1x10-2(31) 

2,9x1 09 

(2,8x1 09) 

5,1  xl  O"2  (3,0x1 01) 

RAF 

C 

1,3x1 011  (5,1x1 09) 

2,5x1 09 

(2,2x1 09) 

2,4  x  1 09  (1 ,5  x  1 09) 

01 

6,2  x  1 06  (1 ,9  x  108) 

2,0x1 09 

(2,0  x1o9) 

6,2  x  1 06  (1 ,8  x  1 08) 

C2 

2,2x  1 08  (3,6x  1 08) 

2,2x1 09 

(2,8  x1o9) 

2,0x1 08  (5,2  x1o1) 

C3 

1,7x10n  (1,9x  1011) 

2,2x1 09 

(2,2  x1o9) 

2,2x  1 09  (2,1  x  1 09) 

C4 

5,9x  1 09  (2,0x  1 08) 

2,2x1 09 

(2,1x1o9) 

1 ,6  x  1 09  (1 ,8  x  1 08) 

C5 

3,4x1 08  (2,4x1 09) 

2,1x1 09 

(2,1  xl  O9) 

2,9  x  1 08  (1,1  x  1 09) 

C6 

1 ,8  x  1 07  (7,0  x  1 08) 

2,2x1 09 

(2,1  x  1 09) 

1 ,8  x  1 07  (5,3  x  1 08) 

SET 

3,1  x  1 09  (9,6x  1 02) 

7,8x1 09 

(8,9  x1o9) 

2,2x  1 09  (9,6x  1 02) 

*  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 
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TS2-H3  2,7  (2,7) 


TS2-H  7,2  (2,7) 


TS2-H5  20,4(17,9) 


TS2-H2  20,6  (18,8) 


TS2-H6  21,1(17,3) 


Figura  S.1 .  Estructura  de  los  TS  que  se  forman  entre  2  y  el  HOO-  en  el  mecanismo  HT  con  el  nivel  M05-2X/6- 
31 1+G(d,p).  Las  distancias  en  Á.  En  la  parte  inferior  de  cada  estructura  se  indica  el  valor  de  la  energía  en  kcal/mol 
con  respecto  a  los  reactivos.  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 
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aci 

avances 

enciéndase 

ingenierías 


TS2-C  4,2  (6,1) 


TS2-C4  6,0  (8,0) 


TS2-C5  7,7  (6,5) 


TS2-C2  8,0  (7,7) 


1,736 


TS2-C1  10,1(8,0) 
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Figura  S.2.  Estructura  de  los  TS  que  se  forman  entre  2  y  el  HOO-  en  el  mecanismo  RAF  con  el  nivel  M05-2X/6- 
31 1+G(d,p).  Las  distancias  en  Á.  En  la  parte  inferior  de  cada  estructura  se  indica  el  valor  de  la  energía  en  kcal/mol 
con  respecto  a  los  reactivos.  Los  valores  entre  paréntesis  corresponden  al  medio  de  etanoato  de  pentilo. 
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Cálculo  de  la  fracción  molar  para  las  especies  del  3,4-Dihidroxibenzaldehído 
en  función  del  pH. 

En  medio  acuoso,  la  fracción  molar  de  las  especies  H2C7H403,  HC7H403"  y  C7H4032~  se 
calculan  por  medio  de  los  valores  de  pKa  según  el  esquema  de  la  Figura  S.3. 


H 


H 


-  H+ 


OH 


pKai 

7.22 


H 


jH 

h2c7h4o3 


OH 


pKa2 

11.40 


cvh4032- 


Figura  S.3.  Disodación  del  3,4-dihidroxibenzaldehído  en  medio  acuoso. 

La  concentración  de  las  especies  de  H2C7H403,  HC7H403_  y  C7H4032_  dependen,  en  este 
caso  del  valor  del  pH  fisiológico  (7.4)  y  se  representan  por  medio  de  las  reacciones  de 
equilibrio  (S.1)  y  (S.2). 

H2C7H403(ac)  H+(ac)  +  HC1H40~(ac);Kal 
HCnHfi;{ac )  H+(ac)  +  C7H4032-(ac);Ka2 

La  constante  de  equilibrio  para  kal  se  representa  en  la  Ec.  (S.3). 


'  ÍH1C1HA] 


(S.2) 


(S.3) 


el  balance  de  masa  da  lugar  a  la  Ec.  (S.4). 

Ka 

resolviendo  para  las  Ec.  (S.3)  y  (S.4)  se  obtiene  la  Ec.  (S.5). 

c»£,n.o,  [  r/"  J  +  Kax 


(S.4) 


(S.5) 
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aci 

avances 

enciéndase 

ingenierías 


donde  el  cociente 


[HfiHfi,] 

-  representa  la  fracción  molar  para  el  H2C7H403 

Ch2c7h4o3 


(X  ),  usando  'as  relaciones:  [H+]=1 0 _pH  y  Ka=1 0 _pKa  para  el  valor  del  pH  fisiológico  se 
tiene  queX  =  0,398.  De  igual  forma,  para  el  HC7H403-  se  tiene  obtiene  la  Ec.  (S.6). 


2  [HC,H,0 3-] 


Ka2 


\[c,HA2'] 

c 

-IV1 

hc7h4o~  L  J 

(S.7) 


de  las  relaciones  anteriores  se  tiene  que  [  HC7H403-  ]  =  CH^H^  -  [  H2C7H403  ]  y 
C  r  u  n  ~  [  HCJ-EC)  - 1  ,  entonces  resolviendo  para  las  Ec.  (S.6)  y  (S.7)  se  obtiene  la  Ec. 

n2c7n4u3  /  4  o  j  1  y 

(S.8). 


1-X 


X 


H2C1Hfil 


K] 


[//*] 


+  Ka « 


(S.8) 


análogamente  se  obtiene  que  XHC7  =0,602.  Para  encontrar  EH2C7H403,  EHC7EH403_  y 
Xr  u  _  se  sabe  que  la  suma  de  las  fracciones  molares  es  igual  a  uno,  de  este  modo  se 

C7H4°32" 

obtiene  que  Xr  H  _  es  1,8  x  10 '5. 

C7H4U32- 
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